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农作物 害虫 预测 预报 的 多 因子 综合 相关 法 


朱 伯 和 承 


(江苏 省 无 锡 气象 站 》 


”摘要 ”农作物 害虫 的 发 生 是 由 各 方面 因子 决定 的 ， 因 此 在 进行 害虫 预测 预报 时 必须 要 考虑 各 方面 尖子 
的 作用 。 
本 文 介绍 多 因 于 综合 相关 的 预报 方法 , 把 预报 量 和 各 预报 因子 按 一 定 标准 化 为 著 干 级 , 用 特征 资料 “0， 
1? 表 示 之 。 经 直接 分 析 各 因子 与 预报 量 的 “ 单 相关 ?建立 预报 方程 。 


引  ” 富 

自然 界 各 种 现象 , 包括 害虫 的 发 生 , 不 是 孤立 存在 的 , 而 是 存在 着 相互 联系 、 相 互 影 
响 、 相 互 制约 的 关系 。 例 如 ,害虫 的 活动 ,与 各 种 生态 因子 (气象 因子 、 地 理 因 子 、 植被 群 
落 、 其 他 因子 等 等 ) 之 间 存 在 一 定 的 联系 。 这 种 联系 如 果 能 用 数学 关系 表示 出 来 , 那么 我 
们 就 可 以 根据 它们 之 则 的 数学 关系 ,由 预报 因子 已 出 现 的 数值 来 估计 预报 对 象 的 数量 ,从 
而 作出 害虫 的 预测 预报 。 

作者 以 前 介绍 的 还 归 估 计 法 ,能 建立 预报 因子 与 预报 量 之 间 的 关系 〈 朱 伯 承 ,1974)。 
但 由 于 人 过 归 估 计 法 需 解 多 元 线性 代数 方程 ， 这 在 预报 因子 众多 的 情况 下 使 用 起 来 不 很 方 
便 。 

这 里 介绍 多 因子 综合 相关 预报 法 。 本 方法 将 预报 量 和 各 预报 因子 按 一 定 标 准 分 为 n 
级 ,用 特征 资料 “0, 1” 表示 。 这 种 方法 最 大 优点 是 无 需 解 多 元 线性 代数 方程 组 ,经 直接 分 
析 各 因子 与 预报 量 的 “ 单 相关 ” 即 可 建立 预报 方程 。 这 是 便于 各 种 害虫 预测 预报 机 构 使 用 
的 方法 之 一 。 


基本 原理 


设 有 产 个 预报 因子 Xs XT2* "Xm 要 预报 的 害虫 要 素 用 y 表示 。 将 MT 2 Xm) 都 分 
为 级, 则 各 因子 与 预报 量 之 间 的 单 相关 由 列 联 表 ( 表 1) 给 出 。 


并 计算 单 因 子 相 关 
pia— G1,2m; hk 1m1, 2 (GD) 
fk. 
将 :个 s 维 单位 行 向 量 
(1 0 0 ... 0) 
(0 1 0 Qa) 
(0 0 0 :.. 1) 


表示 预报 因子 * 与 预报 量 y 的 状态 , 即 : 
当 》 出 现在 1 级 时 ， 嫉 一 1 和 一 0，… 和 一 0， 
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表 1 
OO _ 害 y | 
— -一 一 Ee 
ma ! 2 - | 
1 ni nl 2, “0 ns 歼 1 
2 a ni, 1); 7 ns ea 
区 3 2 了 za 2 ns ns 
5 ns ns ns nss | 7! 
He ni nl, ns a ni 
当 ?y 出 现在 2 级 时 ， 六 一 0， Y2 ] ， 一 0， 
当 ? 出 现在 ?级 时 ， 一 0, 加 一 0， 和 一 
二 Xin 1 xi 一 0 Xis YU, 
xj 出 现在 2 级 时 ,x 一 0, 一] xi 一 0， 
当 x; 出 现在 * 级 时 ， Xi 一 0，xP 一 0， Xi 
第 ;i 个 预报 因子 对 y 一 ! 级 出 现 的 概率 贡献 为 
好 一 2 plexig — plnxa tt phaxiz t+ 十 下 sx (2) 
太 二 1 


对 才 个 预报 因子 ,我 们 取 w 个 因子 概率 贡献 的 平均 值 作为 > 为 1 级 出 现 的 概率 估计 值 , 则 
有 


了 


之 ; 7 DIRXIK (3) 


| 


由 于 各 项 预报 因子 只 能 出 现在 个 级 别 中 的 一 级 ， 因 此 无论 因子 出 现在 哪 一 级 同样 恒 有 
S Xik 1 (4) 
利用 (4) 式 ,可 把 (3) 化 为 


声 二 二 >, Din 十 2, >, 工人 一 PI xin (5) 
?2 ;=1 j=1 k=2 7 
引进 记号 
Di: 一 一 > pi 
1 (6) 
pi 一 本 Cpiig 一 友人 
则 得 预报 方程 
B= ont 2D) or (人 一 1 2 (7) 


了 一 1 =2 
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这 里 要 指出 ，# 是 概率 值 , 必须 满足 
> = 1 (8) 
i=1 


因此 ,至 少 并 仅 需 * 一 1 个 预报 方程 , 第 :级 的 预报 结果 可 由 一 1 一 之 ) 有 得到。 不 过 
i=1 
我 们 不 妨 用 (1)(6)(7) 分 别 建立 个 预报 方程 ,并 由 


Spin=1l KR=1,2s,j=1,2..m) (9) 
i=1 
可 得 
> bn 二 1 | 
(10) 
> 60 
{=1 


也 就 是 说 , (7) 中 :个 预报 方程 的 常数 项 之 和 恒 等 于 1, 同一 变数 项 的 相应 系数 之 和 恒 等 
于 0 ,这 种 关系 可 以 帮助 我 们 检验 预报 方程 建立 得 是 否 正确 。 

预报 : 在 作 预 报时 各 预报 因子 出 现在 哪 一 级 则 取 该 级 为 1 其 他 各 级 为 0, 代入 预报 
方程 算出 各 级 数值 后 , 即 肥 所 值 中 最 大 者 对 应 的 级 别 作 为 预报 。 


预报 举 例 
下 面 举 无 锡 地 区 第 一 代 三 化 晨 蛾 发 生 高 峰 日 期 预报 的 例子 。 历 年 (1960 一 1971 年 ) 
无 锡 县 第 一 代 三 化 螟 峨 发 生 高 峰 日 期 ( 即 预 报 量 y) 由 表 2 给 出 。 
表 2 





年 份 1963 | 1964 196€ | .1967 | 1968 | 1969 | 1970 | 1971 





y 5/23 | 5/20 5/22 | 5/26 | 5/27 | 5/23 | 5/23 | 5/27 





1. 选 择 预报 因子 ”本 例 选 定 如 下 因子 ， 一 一 无 锡 县 1 月 份 雨 量 , 单 位 毫 米 ; x; 
无 锡 县 4 月 份 月 平均 气压 , 单位 毫 巴 ; xs 一 一 无 锡 县 4 月 份 月 平均 温度 , 单位 度 (摄氏 )。 
预报 因子 的 历年 值 由 表 3 给 出 。 











表 3 
















1960 | 1961 





1962 





1963 





1964 | 1965 | 1966 





1969 | 1970 | 1971 


47.5 | 42.9 0.2 | 20.2 | 67.0 74.2 | .15.9 | 26.4 






















14.1 15.4 





13.0 | 14.3 | 16.8 








12.2 14.0 14.5 13.4 13.9 


2. 资料 处 理 : 
(1) 项 报 量 7 分 级 如 下 : 5 月 22 日 以 前 为 1 级 ,5 月 23 一 25 日 为 2 级 ,5 月 26 日 以 
后 为 3 级 。 
《2) 预报 因子 分 级 如 下 ;xi 








当 x 过 10.0 上 毫米 时 化 为 1 级 ,， 当 10.0 二 志 25.0 
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点 米 化 为 2 级 , 当 25.0 瞩 米 过 x 化 为 3 级 。x 当 岂 委 1014.8 毫 巴 化 为 1 级 , 当 1014.8 
一 六 所 1016.0 毫 巴 时 化 为 2 级 , 当 1016.0 毫 巴 二 坟 化 为 3 级 。 ws 一 - 当 吉 三 14.0 度 化 














为 1 级 ， 


闫 报 屋 己 预 报 因 于 的 分 级 值 由 表 4 


1960 | 1961 


“~ 年 份 
因 于 一 


1962 


1963 | 1964 





给 出 。 
表 4 








1965 


当 14.0 之 xs 三 15.0 谋 化 为 2 级 , 当 15.0 度 二 x 化 为 3 级 。 


1966 | 1967 | 1968 | 1969 


3. 分 析 单 因子 相关 ”根据 表 4 求 出 列 联 表 ( 表 5、 表 6.、 表 7)。 


表 5 





p=0, p=0, p=1, pz=0, fiz=1l, p= 0, 
pil 一 了 pils 一 了 p313 一 二 
pih 一 5， pl 一 ， 2 一 0， Phn=0, pb 一 1， 2 一 0， 
pis—= 0， pis = pi 一 - 
pil 一 pe po = pr pal 一 > pin = 二 pal = 二 ， p3l = 0， 
pis—=1, pis=0, ps=0 
4. 建立 预报 方程 ”由 (6) 式 求 出 
pu—=0.3143, ba=0.1143, 6a 一 0.5714 
= 0, bh = 一 0.2667， bi, 一 0.0635 
bis=0.1905,， 6 二 一 0.2667， 号 一 0.2857 
总 一 0.3333， b=0.2667, 5, = 0.1746 
bis 一 0.0476, b= —0.0lll, 2 一 一 4.0476 
2 一 一 0.3333， 5 一 0， bi, 一 一 0.2381 
名 一 一 0.2381， 5353 一 0.2778， 2 = 一 0.2381 





家 6 
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根据 (7), 预 报 方程 为 
fi < 0.3143 十 0.19035x13 0.2667 X22 | 0.26677x23 十 有 十 0.2857x33 


Pp = 0.1143 十 0.3333x1; 十 0.0476x13 十 0.2667x2s 一 0.0111xzs 十 0.1746x3: 一 0.0476xs 
Ps = 0.5714 一 0.3333x 一 0.2381xs 十 0.2778x23 一 0.2381x3 一 0.2381x3s3 (11) 
5. 统计 历史 概括 率 及 预报 ”在 作 预 报时 ,将 因子 按 出 现在 哪 一 级 则 取 该 级 为 1 ,其 他 
级 为 0 ,分 别 代 人 预报 方程 筑 出 记 各 级 数值 , 按 概率 最 大 者 对 应 级 别 作为 预报 。 
根据 所 用 历史 资料 对 12 个 个 例 (1960 一 1971 年 ) 作出 的 预报 列 于 表 8, 与 实 妮 对照， 
12 次 中 报 对 10 次 ,有 2 次 预报 错 1 个 级 别 , 统计 其 历史 概括 率 为 83.3%。 


表 8 































































、 X11 X12 X21 X23 X3t | X32 | X33 Pi. p? pi 预报 | y | 检验 
年 份 八 
1960 0 0 1 1 0 0 0 1 0 10.5683 | 0.3365 | 0.0952 | 1 | 1 1 V 
1961 0 0 1 1 0 0 0 0 1 |0.7905 | 0.1143 | 0.0952 | 1 | 工 1V 
1962 1 0 0 0 0 1 1 0. | 0 10.0476 | 0.1032 1 0.8492 | 3 1 3 1 V 
1963 0 1 0 0 1 0 0 1 0 |0.1111 | 0.8889 0 212 1V 
1964 0 0 1 1 0 0 0 0 1 |0.7905 | 0.1143 | 0.0952 | 1 | 1 1V 
1965 1 0 0 0 0 1 1 0 0 | 0.0476 | 0.1032 1 0.8492 | 3 | 3 | V 
1966 0 0 1 1 0 0 1 0 0 |0.5048 | 0.1619 1 0.3333 | 1 | 1 1Y 
1967 1 0 0 0 0 1 1 0 0 10.0476 | 0.1032 | 0.8492 | 3 | 3 | V 
1968 0 0 1 0 0 1 1 0 0 10.2381 | 0.1508 | 0.6111 | 3 1 3 |V 
1969 0 0 1 1 0 0 0 1 0 |0.5683 | 0.3365 | 0.0952 | 1 | 2 | X 
1970 0 1 0 0 0 1 1 0 0 |0.0476 | 0.4365 | 0.5159 | 3 | 2 | X 
1971 0 0 1 0 0 | 1 1 0 0 |0.2381 | 0.1508 | 0.6L11 | 3 13 |V 


”我 们 用 预报 公式 (I1) 作 了 1972、1973、1974 年 第 一 代 三 化 坚 峨 发 生 高 峰 日 期 的 预报 。 
1972 年 : xi 一 15.3 毫米 , x; 一 1016.0 毫 巴 , x; 一 13.9 度 。 
按 因子 分 级 标准 化 为 : 2z1 一 1, x 一 ry 一 0; xz22 一 1 x 一 x 一 0; wa 一 1, x3; 一 X33 一 
0。 代 入 预报 方程 算出 所 一 0.0476, p 一 0.7143, 所 一 0.2381。 其 中 以 如 为 最 大 ,因此 预 
报 1972 年 第 一 代 三 化 螟 蛾 发 生 高 峰 日 期 为 5 月 23 一 25 日 ,1972 年 实际 出 现 的 日 期 为 5 
月 25 日。 
1973 年 : x1 二 32.9 毫米 ,，z 一 1012.4 毫 巴 ，z 一 15.9 度 。 
按 因 子 分 级 标准 化 为 : X13 一 1， Xl 一 X12 一 0; Yt = 1， X22 = Ya = 0; ra = 1, Ya 一 x 
一 0, 代 人 预报 方程 算出 真一 0.7905, 记 一 0.1143, 包 一 0.0952, 其 中 以 扣 为 最 大 , 因此 
预报 1973 年 第 一 代 三 化 螟 峨 发 生 高 峰 日 期 在 5 月 22 日 以 前 , 1973 年 实际 出 现 的 日 期 为 
5 月 19 日 。 
1974 年 : x 一 57.3 毫米 ,，x; 一 1012.6 毫 巴 , x; 一 15.3 度 。 
按 因 子 分 级 标准 化 为 : x3' = 1, x 一 :xi 一 0， X2 1 Y= 223 一 0; ry 一 1，xat 一 X32 
一 0。 代入 预报 方程 可 算出 名 二 0.7905， 思 一 0.1143， 饭 一 0.0952。 其 中 以 所 为 最 大 ， 
故 预报 1974 年 第 一 代 三 化 晨 蛾 发 生 高 峰 日 期 应 在 5 月 22 日 以 前 , 1974 年 实际 出 现 的 日 
期 为 5 月 19 日 。 
从 1972 一 1974 年 的 预报 看 采 , 三 年 的 级 别 都 预 根 对 了 。 
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讨 论 

据 近 年 来 多 次 使 用 此 方法 作 害 虫 预测 预报 ,我 们 有 如 下 -一 些 体会 。 

1. 本 方法 的 优 缺点 ”优点 是 : 由 于 原始 资料 作 了 (0.1) 化 处 理 , 又 不 需求 解 线 性 代数 
方程 组 ,因此 计算 工作 很 简单 。 由 于 资料 分 为 多 级 。 因 而 预报 的 结果 也 较 细 ,可 以 分 为 多 
种 情 交 作出 预报 。 缺 点 是 : 只 考虑 了 每 个 因子 与 预报 量 的 相关 关系 ， 没 有 考虑 因子 之 间 
父 叉 相关 对 预报 量 的 贡献 。 

2. 关于 挑选 因子 ”由 表 1 可知， 若 某 因子 列 联 表 中 对 角 线 上 的 频数 大 于 其 所 在 行 的 
其 他 各 频数 ,对 于 因子 与 预报 量 正 相关 (预报 举例 中 因子 2) 有 





ne 之 2 或 pin 之 Pi (=D) (12) 
对 于 因子 与 预报 量 之 间 负 相关 (预报 举例 中 的 因子 1、3) 有 
,C1) > nk 或 Pca Dt 一 pilin (ls— 《天 — 1)) (13) 
这 样 的 因子 才 有 预报 意义 。 
男 外 , 阁 预 报 因子 与 预报 量 独 立 , 则 有 
ph 一 PLP (14) 
对 于 正 相关 ,因子 需 符合 
PP > pik (15) 
对 于 负 相 关 ,因子 需 符合 
PDE > pistay) (16) 


究竟 村 大 多 少 ,因子 才 是 可 用 的 9 这 要 视 情 况 及 需要 机 定 。 据 我 们 的 实践 体会 ,一 般 

说 来 20% 左右 就 可 以 了 。 即 正 相关 px 一 pis 之 20% 
负 相 关 pooa-vix 一 站 -ob 人 > 20% 的 因子 是 可 用 的 。 

在 挑选 因子 时 ， 应 尽量 挑选 因子 之 间 相 关 小 的 作为 预报 因子 。 基 两 个 完全 相关 的 因 
于 先入 ,等 于 其 中 一 个 无 用 。 另 外 ,因子 的 个 数 不 宜 太 人 少 , 否 则 失去 综合 考虑 的 意义 ,但 也 
不 宜 太 多 , 太 多 不 仅 增 加 计算 量 , 而 且 会 淹没 主要 因子 的 作用 。 一 般 说 来 3 一 5 个 因子 比 
较 适 宜 。 

3. 关于 定 预 报 临界 值 ”目前 可 用 两 种 方法 : 

(1) 按 最 高 概率 原则 选 定 预 报 临界 值 ( 如 预报 举例 中 所 述 )o 

(2) 以 预报 对 象 出 现 的 平均 概率 为 准 ， 则 可 规定 当 计算 出 的 某 级 概率 大 于 预报 对 象 
间 级 的 历史 平均 概率 时 即 报 该 级 。 

若 预 报 对 象 分 为 多 级 , 当 各 级 个 例 数 相差 不 大 ,而 每 个 级 别 中 有 足够 多 的 个 例 数 使 计 
算 概 率 值 代表 性 较 好 时 ,可 取 概 率 最 大 值 对 应 的 级 别 作 预 报 。 

顺便 指出 ,有 时 也 可 按 预 报 要 求 定 临 界 值 。 如 预报 害虫 大 发 生 时 , 希望 尘 报 要 少 , 在 
定 临 界 值 时 就 可 以 将 预报 害虫 大 发 生出 现 的 临界 值 适 当 取 小 些 。 
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